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W ■ Abstract 

CLc Blazars have been regarded as one of the most powerful sources of the highest 

O ■ energy cosmic rays and also their byproducts, neutrinos. Provided that a magnetized 

'^ ' filamentary system is established in a blazar jet as well, we could apply the mecha- 

nism of multi-stage diffusive shock acceleration to a feasible TeV emitter, Mrk 501 
to evaluate the achievable maximum energy of protons. Taking conceivable energy 
restriction into account systematically, it seems adequate to say that EeV-protons are 
O . produced at this site by our present model. We also estimate neutrino fluxes gener- 

o> ■ 

Q ■ ated by these accelerated protons and discuss the detectability based on an updated 

f^ ■ kilometre-scale telescope such as IceCube. 

Q ' Key words: galaxies: individual (Markarian 501) — magnetic fields — methods: 

numerical — neutrinos: individual (Markarian 501) — shock waves 
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1. Introduction 

The origin of the ultra-high energy cosmic rays (UHECRs) beyond 100 EeV has still been 
an enigma in modern astrophysics. UHECRs are frequently accompanied by significant fluxes 
of energetic neutrinos generated in processes associated with hadronic interactions. Although 
the source direction of cosmic rays with energy above 60 EeV can be identified by an analysis 
of recent Auger data (Abraham et al. 2007), neutrinos are regarded as being the most useful 
probe of their source confirmation owing to no deflection by galactic/intergalactic magnetic 
fields and the weakest reduction by interactions with the cosmic background radiation. The 
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models for the origin of the highest energy cosmic rays are classified into two categories: top- 
down and bottom-up. In top-down scenarios, neutrinos are presumed to be decay or annihilation 
products of cosmological remnants with the Grand Unified energy scale of ~ 10^^ GeV. Within 
this framework, numerous models have been proposed so far, such as topological defects, Z- 
bursts, annihilation of dark matter particles and so on (for a review, see Bhattacharjee & Sigl 
2000). In an alternative approach, the bottom-up model is based on the idea that the highest 
energy cosmic rays are generated by astrophysical objects. The most plausible candidates are 
gamma ray bursts (GRB: Waxman & Bahcall 1997; Vietri 1998) and active galactic nuclei 
(AGN: Honda & Honda 2004a; Honda 2009). In addition to these objects, it is pointed out 
that UHECRs are also produced via stochastic acceleration in the giant lobes of radio galaxy 
(Cen A: Fraschetti & Melia 2008; O'Sullivan et al. 2009), which accounts for the past year 
discovery by the HiRes and Auger collaborations. Supernova remnants. X-ray binaries, mini- 
quasars (e.g. Gaisser, Halzen, & Stanev 1995, for a review), and soft gamma ray repeaters 
(SGR: Halzen et al. 2005; loka et al. 2005) can also be stellar type neutrino sources, any of 
which could produce observable fluxes of energetic neutrinos. If protons are accelerated at these 
astrophysical objects, neutrinos are expected to be produced in collisions with ambient photons 
or protons. 

Some of the above-mentioned models predict neutrino flux at the level of a few events per 
km^ year. The upper range of this estimation seems to be within reach of a first-generation neu- 
trino telescope, such as AMANDA (Antarctic Muon And Neutrinos Detector Array). Among 
a wide variety of physics topics to be explored with the neutrino telescope, the most important 
goal is to search for the origin of cosmic radiation, especially that originated from AGN and/or 
GRB. For this purpose, AMANDA has been primarily optimised in the energy range from TeV 
to PeV, targeting not only the diffuse, but point source, flux of energetic neutrinos. According 
to an analysis of AMANDA data, we unfortunately have no evidence for point sources so far, 
(Woschnagg et al. 2005; Abbasi et al. 2009), except for a temporal coincidence with an or- 
phan flare of the TeV blazar, lES 1959+650 (Halzen & Hooper 2005). The IceCube, a cubic 
kilometre-scale neutrino telescope, which is partly operating and still under construction as a 
successor to AMANDA at the same site, is designed to detect the fluxes about 10-50 events 
per km^ year. Similar event rates are predicted by model calculations assuming that AGN or 
GRB are the actual source (Halzen 2005) and hence the IceCube is expected to certainly detect 
neutrino fluxes from such point sources. In order to observe a guaranteed source, however, even 
the km^-sized detector should safely operate over a period of at least ten years. 

As for the cosmic accelerators, the most promising mechanism to lead highest energies 
and power-law energy spectrum is considered to be of diffusive shock acceleration (DSA: Drury 
1983; Lagage & Cesarsky 1983a; 1983b). Solar energetic particles and galactic cosmic rays 
are consistently explained by this mechanism applied to shocks at heliosphere and supernova 
remnants, respectively (for a review, see Blandford & Eichler 1987); these acceleration sites have 



actually been confirmed by some observations (Koyama et al. 1995). Polarization data in the 
range of radio to optical wavelengths imply that in situ acceleration of electrons (and possibly 
ions) must also be taking place in the knots of AGN jets (Honda & Honda 2004a) or in the hot 
spots of Fanaroff- Riley type-II (FR-II) sources (probably terminal shocks) (Perley et al. 1984; 
Carilli & Barthel 1996). At these sites, magnetic fields in the vicinity of a shock front plays an 
essential role in particle acceleration. A detailed configuration of the fields in extragalactic jets 
has been revealed by polarization measurements using very long-baseline interferometry. For 
example, the quite smooth fields predominantly transverse to the jet axis, are typically observed 
in the core region of BL Lac objects (1803+784: Gabuzda 1999; 0300+470: Nan et al. 1999). 
Allowing the propagation of a shock wave through the jets, this implies the establishment of 
a (quasi-) perpendicular shock, which is capable of accelerating cosmic ray particles efficiently. 
Evidence for large-scale toroidal magnetic fields has been also discovered in the galactic center 
(GC) region (Novak et al. 2003). Based on the above results, it seems adequate that system of 
magnetized filaments has been established in blazar jets as well. 

We have suggested a theoretical model to account for the generation of such a large-scale 
toroidal (transverse) magnetic field in astrophysical jets (Honda & Honda 2002). In this model, 
huge currents launched from a central engine are broken into many filaments whose transverse 
sizes are self-adjusted for the electromagnetic current filamentation instability (CFI: Honda 
2004 and references therein). In the nonlinear stage of the CFI, the magnetized filaments are 
often regarded as being strong turbulence that can largely defiect charged particles. Allowing 
shock propagation through the jet, the particles are expected to be quite efficiently accelerated 
for the DSA scenario, which appears to be favorably taking place in the AGN jets. Indeed, some 
knots of a radio galaxy jet are associated with the shock fronts (M87: Biretta et al. 1983; Capetti 
et al. 1997) and circumstantial evidence for in-situ acceleration of electrons (Meisenheimer et 
al. 1989) have been found. It was also pointed that a similar pattern of small-scale quasi- 
static magnetic fields can also be established by some numerical simulations: e.g., during the 
collision of electron-positron plasmas existing in SNRs, pulsar winds, GRBs, relativistic jets, 
and so on (Kazimura et al. 1998; Silva et al. 2003). Using a three-dimensional relativistic 
electromagnetic code, Nishikawa (2003) show that non-uniform small-scale magnetic field is 
generated due to the Weibel instability at a jet front propagating through an ambient plasma 
with/without initial fields. These dynamics might be involved in some of the knots in the FR-I 
radio jets, which appear to be a shock established when faster material is overtaking with slower 
one (M87: Biretta et al. 1983). 

In the case of blazars, however, the detailed internal structure of their jets has still not 
been confirmed because of optical thickness with respect to the observer's line of sight. Another 
approach to shed light on their configuration is provided by the remarkable short variability 
timescale of a blazar light curve (e.g., Mrk 421: reaching a few minutes, Cui 2004; Blazejowski 
et al. 2005). This is significantly shorter than the light-crossing time at the black hole horizon. 



which iinphes the presence of some substructure in the parsec-scale jet. According to our 
model configuration of the filamentary jet, the strength of the local magnetic field is described 
by a power-law of the filaments transverse sizes. If charged particles are injected into this 
system, they are diffusively accelerated by a collisionless shock being scattered by the field 
fiuctuations. Since the efficiency of the acceleration and loss depends upon the spatial size 
scales, the local maximum energies of accelerated particles are also characterized by their nearest 
filament sizes. The spectrum extending to the X-ray region is attributed to the synchrotron 
radiation of accelerated electrons. In particular, the correlation between X-ray and gamma-ray 
light curves of Mrk 421 is well reproduced by a model of the hierachical turbulent structure of 
the jet (Honda 2008). The most interesting consequence is that the transition of a hierachical 
turbulent structure seems to be responsible for the observed patterns of energy-dependent light 
curves, e.g., soft lag (Takahashi et al. 1996; Rebillot et al. 2006) or hard lag (Fossati et al. 
2000) and a tight-correlation mode (Senibay et al. 2002). 

As for Mrk 501, which is also firmly established as a TeV 7-ray emitter, the X-ray light 
curve shows a very rapid fiare varying over several 100 seconds (Catanese & Sambruna 2000). 
It is also revealed that the TeV fiares correlate with X-ray radiation on timescales of hours or 
less from multiwavelengths observations (Plan et al. 1998; Catanese et al. 1997; Krawczynski 
et al. 2000; Sambruna et al. 2000). Moreover, the largest shift of the peak energy during the 
peak-luminosity change was observed among all blazars (Kataoka et al. 2001). Summarizing 
these observational results, we can say that the acceleration of electrons similar to Mrk 421 is 
also taking place in Mrk 501. Presuming a finite ion abundance in the jet of Mrk 501 (Rawlings 
& Saunders 1991), such a DSA mechanism operates for arbitrary nuclei (of course including 
protons) as well. Indeed, from the X-ray spectrum of SS 433 jet, emission lines of various 
elements such as Ne, Mg, Si, and Fe were observed (Kotani et al. 1996). 

In the present paper, we evaluate the achievable maximum energy of protons in the Mrk 
501 jet to estimate the resultant neutrino fiuxes in accordance with our diffusion and accelera- 
tion mechanism in a filamentary jet, which we have proposed and developed in a series of papers 
(Honda & Honda 2004b; 2005; 2007). In section 2, installing our model structure of magnetized 
filaments to the jet of Mrk 501, we describe the acceleration proccesses of protons using our 
novel DSA mechanism. Taking competitive energy losses and restrictions into consideration, 
we present the scaling of the maximum energy. In §2.1 we describe the acceleration taking 
place at local magnetic fields induced by the current filament each (referred to as preliminary 
acceleration) and in §2.2 the additional acceleration due to the inter-filaments' defiection. In 
§2.3, we calculate the achievable maximum energies of accelerated protons with respect to the 
transverse filament sizes and magnetic field parameters. We then evaluate neutrino fiux from 
Mrk 501 and discuss its detectability based on a km-scale neutrino telescope such as IceCube 
in section 3. Conclusions are summarized in section 4. 



2. Proton Acceleration in Magnetized Filaments 

Suppose a parsec-scale blazar jet transporting energetic particles from the central engine. 
The directional plasma flow will favorably induce huge currents (Appl & Camenzind 1992) 
that breaks up into many fllaments because of electromagnetic CFl. Each fllament generates 
a small-scale transverse magnetic fleld whose strength is determined by the transverse size of 
a filament, A, i.e., \B\ = B^{X/DY'^~'^^/'^, where B^ = \B\x=d; D and /3 = 4.3 denote the jet 
diameter and power-law spectral index of the magnetized filamentary turbulence, respectively 
(Honda & Honda 2007). Viewing the jet globally from outside, randomly oriented fields inside 
jet are cancelled, except for a large-scale toroidal magnetic field along the envelope. 

As was illustrated in our previous paper (Honda & Honda 2005), there exists an energy 
hierarchy for protons trapped in such a filamentary current system: (i) E^ <^ |eA| and (ii) 
Ep ^ |eA|, where e and A are the charge of an electron and the vector potential, respectively. 
The former corresponds to a low-energy regime, in which protons are strictly bounded for the 
local magnetic field induced by each filament. The latter refiects the higher energy regime in 
which the validity condition of the quasilinear approximation, (/p) ^ |5/p| is satisfied, where 
(/p) and 5/p are the averaged and fiuctuated part of the momentum distribution function for 
protons in the test-particle approximation, /p, respectively. In this regime, protons are no 
longer trapped with local magnetic fields, but are deviated with small fluctuations. Below, we 
describe the acceleration and competitive energy losses in each energy hierarchy. 

2.1. Preliminary Acceleration 

In the low-energy regime of E^ <^ |eA|, protons are strictly bounded and gyrating around 
the local magnetic fleld induced by each fllament. Energies of injected protons are elevated via 
the conventional diffusive shock acceleration being resonantly scattered from small magnetic 
fluctuations (Drury 1983). The characteristic acceleration time of protons is described by 

Ve.acc — (^ ^— j \r~r~) ' ^ 

where ?7p = (3/26)(A/2rg_p)^/^, presuming Kolmogorov turbulence, h {<ti 1) is the energy density 
ratio of fluctuating to local mean magnetic flelds, rgp is the proton gyroradius, and r (=4 
for the strong shock limit) is the shock compression ratio. The achievable maximum energies 
of protons are limited by the various radiative cooling timescales, and here we consider two 
representatives: the synchrotron radiation losses, tp,syn, and collisions with ambient photons, 
tp^. In blazar jets the propagation time of shock, tgh, also restricts the highest energy, and 
hence the temporal limit is given by 

''pre, ace I^lll(^^p7; ^p,syn) ''shj- \^) 

Besides the temporal limit, there exists a spatial limit: the maximum gyroradii of protons are 
not allowed to exceed the transverse sizes of in-situ fllaments. The achievable energy of protons 



via the preliminary acceleration is determined by comparing the temporal and spatial limit. In 
the following, conceivable energy constraints for Mrk 501 are provided. 

2.1.1. Synchrotron cooling losses 

Gyrating around the filament- induced local fields, and being scattered by small fluc- 
tuations, protons emit synchrotron photons, which can be a dominant cooling process when 
the energy density of the magnetic field is sufficiently greater than that of the radiation. The 
characteristic loss time for proton synchrotron is given by (Rybicki & Lightman 1979) 

where nip and Trie are the masses of a proton and of an electron, respectively; o"t is the cross 
section for Thomson scattering and 7p is the Lorentz factor of accelerated protons. The average 
energy density of the local magnetic field is denoted by mb = -B^/(87r), where B =\B\. It should 
be noted that the adopted field value of B^ = 0.02 G (Tavecchio & Maraschi 2001) is not the one 
averaged over the compact blob but the maximum for a filament whose radial size is compared 
to the jet diameter. 

2.1.2. Collisions with ambient photons 

We consider proton-photon collisions leading to a pion-production cascade. Since the 
collision timescale is characterized by the target photon spectrum, which is still unknown, we 
adopt here a description by a single power-law: n{e^) oc e^"^ (Bezler et al. 1984), where n{e^) 
is the number density of photons per unit energy interval de^. Then, the timescale can be 
expressed as 

tp^ = X'^ f^j tp,syn ^ 6.0 X 10^° f — j S, (4) 

where m^ ~ 4.0 x 10~^ ergcm"^ (Bicknell et al. 2001) is the average energy density of target 
photon fields and x — 200 (Biermann & Strittmatter 1987). Since u-b/u^ — 10~^, it is found 

XllciL ''p'7 *^^^^ ''p.syn* 

We should also check whether collisions with particles (especially, with protons) becomes 
effective or not. The characteristic cooling time of relativistic protons due to inelastic pp 
collisions can be written as (Aharonian 2004) 

1 ^ ^ .^15 /I cmr^\ /40mb\ 

noappfc \ no J \ (Tpp J 

where the number density of the target hydrogen medium, rig — 1 cm~^, corresponds to the 
upper limit assuming mass loaded hadronic jet models. The averaged total cross section, o"pp, 
is approximately 40 mb at very high energies, and / ^ 0.5 is the coefficient of inelasticity. It 
would be fair to say that collisions with protons should not be taken into consideration. 

2.1.3. Propagation time of shock 

The achievable maximum energy of protons is also restricted by the duration that the 
shock is propagating through the jet from the central engine to the working surface. This 
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corresponds to the age of a blob, which is propagating with mild-relativistic speed through the 
relativistic jet (Honda & Honda 2004a), and is approximately written by L/U . For Mrk 501, 
we obtain 



L ^ . _o/ L \ /0.5c 



where L ~ 7.4 pc is the distance from the core to the blob (Giovannini et al. 1999) and U ^ 0.5c 
is the average speed of the blob (Miicke & Protheroe 2001). This blob is currently operating 
for particle acceleration, and therefore tgh cannot be compared to the adiabatic loss timescale 
for AGN, which corresponds to the lifetime of the shock (Miicke et al. 2003). The timescale of 
adiabatic losses is expressed as tad = 3L/(2rjfi.), where Fj and v^ are the Lorentz factor of the 
jet and the speed of the radial expansion. Considering narrow opening angles of blazars, which 
implies v^ <^ f/, we can safely say that tgh is sufficiently shorter than tad- 
2.1.4- Escape from local filaments 

Once accelerated particles acquire sufficient energies, they are escaping from the local 
magnetic fields, given as a function of the filament sizes. This spatial limit is given by the 
assumption that the gyroradius of each particle should not be beyond the transverse size of the 
local filament. Then, the local maximum energy is denoted by E* = ei?rg p < eBX/2. 

The transverse sizes of the minimum and maximum filaments are regarded as being 
comparable to the Debye sheath and the jet diameter, respectively. Thus, the energy of particles 
escaping from the largest filament can be written as 

P 2-2 V0.02Gy V0.067pcy ' ^ ^ 

where D ~ 0.067 pc is the extent of the radio emitting region (Giroletti et al. 2008). This value 
is estimated by substituting 0.02 G for B in the formula suggested by Marscher (1987). 

2.2. Additional Acceleration 

Once protons are accelerated and injected into the high-energy regime of Ep ^ \eA\, they 
are no longer bounded to the local magnetic field. Field vectors are randomly oriented in the 
internal jet in the transverse direction to the filament, and hence these protons can move almost 
freely in the jet being defiected by the random magnetic field. Allowing for shock propagation, 
protons are additionally accelerated in this regime, being off-resonantly scattered in a forest 
of magnetized filaments. In this aspect, the preliminary accelerator via the conventional DSA 
mechanism can be regarded as being an 'injector' to this further booster. The new injection 
mechanism is referred to as transition injection, in analogy to bound-free transition in atomic 
excitation. The characteristic time of the additional acceleration for protons is given by 

, ^ 3v/67r/3r(r + l) cE^ 

tadd,acc g^^ _ ^^^^^ ^ ^^^^^ ^ 2)(r - 1) e^B'^sDU^ ' ^""^ 

where /3p = 3 (e.g., Stecker & Salamon 1999; de Marco, Blasi, & Olinto 2003) and Ep = 7p?7ipC^ 



are the spectral index and energy of protons, respectively; /3 is the filamentary turbulent spectral 
index, and |i?cffP/-B^ ~ 0.5 is assumed. The maximum energy boosted by the additional 
acceleration is also restricted by some radiative loss processes. In addition to two possible 
timescales: for diffuse synchrotron, td,syn, and for the p7 collision, td,p7, we consider the shock 
propagation timescale, tsh, and therefore tadd,acc = niin(td,syn,^d,p7,4h)- The spatial limit will 
also be considered. 

2.2.1. Diffuse synchrotron losses 

The additionally accelerated protons deflected by the random magnetic fields emit un- 
polarized synchrotron radiation. Here, we adopt the characteristic cooling time for the diffuse 
synchrotron from our previous paper (Honda & Honda 2007), which was derived from the 
theoretical basis by Toptygin & Fleishman (1987): 

3mpc fm^Y 1 (3{(3 + 2Y{(3 + 3) _^ 
'''^^""4aTMB7pVmey SOvr^ 2/5(/32 + 7/3 + 8) " 'p'^^-"^^^^' ^^^ 

where tp^syn is the timescale of normal synchrotron losses and f(/3) = (67r)~^[/3(/3 + 2)^(/3 + 
3)/2^{P'^ + 7/3 + 8)]. For /3 = 4.3 (Honda & Honda 2007), the maximum energy of a proton, 7p, 
is directly estimated by the equation tadd.acc = ^d.syn, which gives 

7p ~ 2.7 X 10" (-^] ' ' ( ^ ] . (10) 

'^ \0.5cJ V0-067pc; ^ ^ 

It should be noted that the maximum energy restricted by the diffuse synchrotron is free from 

the magnetic field strength because of the same S-dependence of tadd.acc and td,syn- 

2.2.2. Energy losses due to jry collision 

Cooling due to collisions with photons is also characterized by the above-mentioned 
diffuse synchrotron: td,p7 = X~^{'^B/u-y)td,syn- Using the same target photon spectrum and 
other notation with those given in §2.1.2, we obtain 

.„. 4.6X10- (5^) (jj^) . (11) 

2.2.3. Restriction from the Age of the Blob 

As we have mentioned in §2.1.3, the shock propagation time regarded as the age of the 
blob restricts the acceleration. Equating tadd,acc = ^h, we obtain 

,„.4.3xlO»fiL)"7^)"7^^)"'^ (12) 

"^ \0.5cJ \7.^pcJ \^0.067pc/ ^ ^ 

2.2.4. Spatial Limit for Additional Acceleration 

It should be noted that additionally accelerated particles are not gyrating to the local 
fields, but deflected by randomly oriented fields. Thus, the spatial limit is determined by the 
condition that the mean free path for this coUisionless deflection should be smaller than the 
system size: £ ~ pD, where p > 1 is a dimensionless parameter (see, Honda & Honda 2005, for 
the details). It gives 
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2(/3-l)(/3p + 2)(/3p + 4)^f ^/ 
\ n P ei^cff-L^ 



as / ^cflf \ I D 



TTTnax 
^P = 

-2.1xl0^«( ^"" )[ 1 eV, (13) 

where p = 1 is assumed. 

2.3. Achievable Maximum Energy of Protons 

Putting together the above discussion, provided botii the normal and transition injection 
of protons are reahzed in a blazar jet, multi-stage diffusive shock acceleration can take place. 
In this subsection, we present the results of a numerical evaluation of the achievable maximum 
energies of protons in both stages within the conceivable parameter region. In Fig.l, the 
energies of accelerated protons are plotted against the transverse filament size parametrizing 
the ratio of the mean/fluctuate magnetic energy density, h. The transverse filament sizes. A, 
covers the range from 10""* pc to 0.067 pc. The minimum and the maximum sizes are compared 
to the correlation length of Alfvenic fluctuation and the jet diameter, respectively. As shown, 
in the case of preliminary acceleration the achievable maximum energies of protons are mainly 
restricted by the escape from individual filaments, which is represented by a linearly increasing 
line (oc A^'^^). In the case of 6 = 10"^, the proton energies achieved via preliminary acceleration, 
E* are determined by r^^ < A over the whole range of the considered filament sizes. In the 
case of 6 = 10^'^, the values of E* are also determined by r^^ < A for A < 1.7 x 10^^ pc, while the 
£'p-values are more severely restricted by the propagation time of the shock for A > 1.7 x 10^^ 
pc, and hence the achievable energy via the preliminary acceleration, E* is cut-off at 10^^ eV. If 
h < 10~^, the -Ep-values are perfectly overlapped with those of 6 = 10~^. On the other hand, the 
maximum energy achieved by the additional acceleration E^'^^ is independent of not only the 
filament size, but also the 6-parameter. The energies of accelerated protons can be additionally 
boosted up to the value determined by the spatial limit expressed with a horizontal line of 
^max = 2.1 X 10^*^ eV, irrespective of their energies acquired by the preliminary acceleration. 

We also present the A-dependence of -E* parametrizing the -Bm-value in Fig. 2. The lines 
corresponding to rg^p = A shift upwards with a larger B^. For the cases of i?m = 0.01 G and 0.1 
G, the proton energies achieved by the preliminary acceleration, E* are determined by rgp < A 
over the whole range of conceivable filament sizes. In the cases of B^ = 1 G and 10 G, the 
protons dominantly lose their energies via collisions with photons for A > 4.2 x 10~^ pc and for 
A > 1.5 X 10~^ pc, respectively. Thus, the resultant maximum value of E* is nearly 10^^ eV. 
As for the additional acceleration, the £'™'"-values are expressed by horizontal lines because of 
their independence of A. The lowest and the highest E™^^ correspond to those of B^ = 10~^ G 
and Sm = 10~^ G, respectively, and in both cases the maximum energies are restricted by the 
spatial limit. The ii^^'^-value for i?^ = 1 G is identical to that for B^ = 10 G, and in these 
cases the energy is restricted by collisions with ambient photons. 

9 



10" 



en 
u 
i> 

e 
o 

o 
10^ 



10' 







A 


Additional Acceleration: E. """' >/ 

P yT 








y/^ b =10"' 




y^ 




y^ 




yT Preliminary Acceleration 




/^=i»' /) = 0.067 pc^ 




1 





JQ-4 JO-3 JQ-2 JQ-1 

Transverse Filament Size [pc] 



Fig. 1. Maximum energies of accelerated protons versus the transverse filament sizes for various 6- values. 
The value of i?„i is fixed at 0.02 G. The monotonically increasing line and its branch denote the maximum 
energies achieved by the preliminary acceleration. Solid and dashed lines correspond to the cases of 
h = 10~^ and lO"'^, respectively. The horizontal bold line denotes the maximum energy achieved by 
additional acceleration. 
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Fig. 2. Maximum proton energies versus the transverse filament sizes for various _Bni- values. The value 
of h is fixed at fO~'^. Four monotonically increasing lines denote the limitations obtained by preliminary 
acceleration. The cases of B^ — 0.01 G, 0.1 G, 1 G, and 10 G are plotted by solid, dotted, dot-dashed, 
and dashed curves, respectively. Horizontal thick lines represent the cutoff energies achieved by additional 
acceleration. The line types are the same as those for preliminary acceleration. 
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3. Estimation of the Neutrino Flux from Mrk 501 

In the previous section we calculated the achievable maximum energies of protons in 
a blazar jet. Since the proton energies accelerated by multi-stage DSA are sufficiently higher 
than the neutrino-producing threshold via photopionization, we can safely expect significant 
fluxes of energetic neutrinos. 

3.1. Neutrino Energies Produced by Accelerated Protons 

In the photomeson process of p7 — t- A^^ -^ n^N, the center-of-mass energies of acceler- 
ated protons should be beyond the A-resonance threshold, m^ = 1230 MeV, 

2 2 

^mm _ p2 A p .^^ 

where Fj and E^ are the Lorentz factor of a blazar jet and the mean energy of ambient photons, 
respectively. Presuming that the generated pion energy is equally divided into four leptons 
through the decay 7r+ — )■ u^fi^ — t- u^e'^i'f.Uf,, the neutrino energy is described as 

E^ = -{xp^^)Ep, (15) 

where (xp^T^) ~ 0.2 is the averaged fraction of momentum transfer from a proton to a pion. 

As for the Mrk 501, the peak energy of synchrotron emssion in the X-ray band is shifted 
from 1 keV to 100 keV during flaring (Kataoka et al. 2001). Since these energies are achievable 
for synchrotron photons emitted from co-accelerated electrons, we take these values as the mean 
energy of ambient photons. For E^ = 100 keV, we obtain the neutrino producing threshold using 
eq.(14) as 

TiXVlOOkeV^ 




Ey 



PeV, (16) 



which leads to Ef"" = 70 TeV. Similarly, for E^ = 1 keV, Ef"" = 140 PeV and hence Ef"" = 7 
PeV. Since the maximum energy of accelerated protons is sufficiently higher than these threshold 
values, Mrk 501 can be a feasible source of neutrinos. In both cases we determined E'^'^^ ~ 110 
PeV, according to the upper limit estimated in §2.3. 

3.2. Energetic Neutrino Flux from Mrk 501 

In order to calculate the neutrino flux produced at an individual blazar, we begin with 
the following useful formula: 

'~E^^dE^ = ^, (17) 

£"'" dE^, Andi 

where the upper and lower limits of integration are the maximum and minimum (threshold) 
neutrino energies, estimated in §3.1, respectively. The luminosity distance to the blazar, di, is 
deflned as 



11 



rfL = <(i + z) = "S^ilA r J(i + zy{i + n^z)-z{2 + z)nAdz, (is) 

-no "'0 ^ 
where z is the redshift (=0.034 for Mrk 501: Quinn et aL 1996) for Mrk 501 and d^ is the 
proper motion distance, which is written by the Hubble constant, Hq = 71 km~^Mpc~^, the 
matter density, Q^i = 0.27 ±0.04, and dark energy density, Q^ = 0.73 ± 0.04, adopted from 
WMAP results (Bennett et al. 2003). 

The observed neutrino luminosity, L^, is defined as 

L.^^, (19) 

where A^;^, (E^), and Atobs — Id are the number of produced neutrinos, the mean neutrino 
energy, and the observed variability time of the emission region, respectively. Since N^ and 
(El,) are not observables, we introduce the optical depth defined as 

r = Rn^ajry^^ = R CTp^^A, (20) 

V{Ey} 

where L^ ~ 10^^ erg/s is the observed photon luminosity, V represents the volume of emission 
region with radius R, (E^) ~ 10 eV denotes the mean energy of photons, and o"p^->A — 10~^* cm^ 
is the A-resonance photopionization cross section. Then, we can express L^ with observable 
quantities in a compact notation, 

L, = KrLp = Kre^Lp^ob,, (21) 

where K ~ 0.024 is a constant taking the branching ratios of interaction chain and corresponding 
momentum transfer into account. The intrinsic proton luminosity, Lp, can be replaced with 
the observed Lp^obs taking the proton interactions with ambient photons into consideration. 
However, Lp^obs is still unknown; we assume it to be as 10 % of the photon total luminosity 
(Halzen & Hooper 2002). We also replace r with (1 — e""^), with considering the possible p and 
TT"*" absorption in the vicinity of the source. 

Summarizing the above discussion, we adopt the following formula of differential flux 
from an individual blazar (Bazo & Gago 2005): 

o 0.1A'L-^(l-e-^)e(i-'=~^) 

dK j^!::E.E;;^e.p(--^dE^'"^ '"H E^^V' ^ ^ 

where the factor 2 reflects the contribution of neutrinos and anti-neutrinos. Since we are 
concerned with the point-source neutrino flux, we consider the i/^ contribution on account of 
their highest angular resolution. This is attributed to the large muon track, having advantege 
of reconstructing paths. 

In Fig.3 we plot the differential flux of neutrinos from Mrk 501 against their energies. The 
energy of neutrinos is limited to the range between the maximum and the minimum predicted 
by this model. Here, we adopt the typical power-law energy distribution of neutrinos, E~p, 
with the exponential cutoff, E^^ . We parametrize the power-law index p = 1.3, 2.0, 2.7, where 
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Fig. 3. Differential neutrino flux from a blazar jet of Mrk 501 versus neutrino energy. The power-law 
energy distribution of neutrinos, E~p, with an exponential cutoff E^^ = 2.1 x 10"'^^ eV is assumed. The 
power-law indices of p = 1.3, 2.0, 2.7 are denoted by the dashed, dotted, and solid curves, respectively. 

p = 1.3 and p = 2.7 are taken from the upper and lower limits of the exponents, respectively, and 
p = 2.0 is a resemblance of the parent proton spectrum supposed to be accelerated by Fermi's 
mechanism. The cutoff energy, E^^ = 2.1 x 10^ GeV, is the achievable maximum energy of 
AGN protons estimated in §2.3. One can see that for a steeper energy spectrum the neutrino 
flux is more significantly reduced in the energy range of E^, ^ 10^ GeV. On the contrary, in 
the lower energy region of E^ ^ 10^ GeV the flux is higher in the case of a steeper spectrum 
because of the contribution from the denominator. Anyhow, the neutrino flux originating from 
a single blazar is trivial compared to the observable level. 

3.3. Calculation of the v ^-Induced Event Number 

It is more difficult to detect higher energy neutrinos because of the steeply falling spec- 
trum. In order to detect the interaction of a TeV neutrino with a Cherenkov telescope, the 
effective volume of the detector is required for kilometer scale to cover the typical range of TeV 
muons. In the effective region of the telescope, the probability to detect z/^ in the TeV-PeV 
range is approximately given by 



P. 



R,, 



1.3 X 10" 



E. \ 



0.8 



(23) 



Aint VlTeVy ' 

where R^ and Amt are the muon range and the neutrino interaction length, respectively (Gaisser, 
Halzen, & Stanev 1995). 

We can now compute the diffuse neutrino event rate by integrating the total differential 
flux multiplied by the detection probability of eq.(23) 
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_£;min 



dE,, 



-{k;, 



p, 



Ufi-i-tii 



(24) 



tot 



where the upper and lower hmits of the integral are E^^^ ^ 1.1 x 10^^ eV and ii^™™ ^ 7 x 10^^ eV 
calculated in § 3.1, respectively. In order to compute the diffuse neutrino flux from the observed 
blazar distribution, we should integrate the flux from all blazars while taking the Doppler factor 
distribution into consideration. For simplicity, here we adopt the effective number of blazars 
(Halzen & Zas 1997), using the 7-ray flux ratio of diffusive to the single blazar, which is derived 
from the luminosity function of 20 brightest blazars obtained by the EGRET collaboration 
(Chiang et al. 1995). Thus, instead of summing up the total isotropic differential flux of 
neutrinos ( ^^^^ ) , we simply estimate the diffuse flux by multiplying the flux of Mrk 501 

\ "'" / tot 

by the resultant number 130 sr"^. A correction for the difference in the spectral indices of 
gamma-ray and neutrino fluxes enhances the neutrino flux by a factor of three. We present the 
estimated neutrino event rates in Table 1 for conceivable power-law exponents. Considering 
the 4:71 coverage and typical exposure time (~ 10 yr), the value N^ ^ 2 km~^ yr~^ sufficiently 
reaches to the observable level. 



Table 1. Neutrino event rates for various power-law exponents. 



Power-law exponent 


Event rate 


p=1.3 


0.78 


p = 2.Q 


1.3 


p = 2.7 


1.9 



4. Conclusions 

In this paper we applied our model of multi-stage diffusive shock acceleration to the 
blazar Mrk 501 and evaluated the maximum energies of protons. We also calculated the neutrino 
flux from this single source generated by these accelerated protons, and discuss the detectability 
with the updated neutrino telescope, IceCube. We obtained the EeV protons and differential 
flux during the operating time of the detector by not a simple, optimistic order estimation, but 
by considering the systematical energy restriction at the source. The results are rather sensitive 
to the values of the magnetic field strength and the size of the accelerator. The S-field value is 
related to the system size via the radiation intensity of the emission region. In the conceivable 
range of the field strength, EeV protons are obtained by additional acceleration. 

Y. S. H. thanks K. Murase for useful comments. 



14 



References 

Abbasi, R. et al. 2009, Phys. Rev. D, 79, 062001 

Abraham, J. et al. (The Pierre Auger Cohaboration) 2007, Science, 318, 938 

Aharonian, F. A. 2004, in Very High Energy Cosmic Gamma Radiation: A Crucial Window on the 

Extreme Universe (Singapore: World Scientific), 109 
Appl, S., & Camenzind, M. 1992, A&A, 256, 354 
Bazo, J. L., & Gago, A. M. 2005, arXiv: astro-ph/0504554 
Bennett, C. L., et al. 2003, ApJS, 148, 97 
Bezler, M., Kendziorra, E., Staubert, R., Hasinger, G., Pietsch, W., Reppin, C, Triimper, J., & 

Voges, W. 1984, A&A, 136, 351 
Bhattacharjee, P., & Sigl, G. 2000, Phys. Rep., 327, 109 

Bickneh, G. V., Wagner, S. J., & Groves, B. 2001, in Proc. AIP Conf., 558, 261 
Biermann, P. L., & Strittmatter, P. A. 1987, ApJ, 322, 643 
Biretta, J. A., Owen, F. N., & Hardee, P. E. 1983, ApJ, 274, L27 
Blandford, R., & Eichler, D. 1987, Phys. Rep., 154, 1 
Blazejowski, M., et al. 2005, ApJ, 630, 130 

Capetti, A., Macchetto, F. D., Sparks, W. B., & Biretta, J. A. 1997, A&A, 317, 637 
Carilh, C. L., & Barthel, P. D. 1996, A&A Rev., 7, 1 
Catanese, M., et al. 1997, ApJ, 487, L143 
Catanese, M., & Sambruna, R. M. 2000, ApJ, 534, L39 

Chiang, J., Fichtel, C. E., von Montigny, C, Nolan, P. L., & Petrosian, V. 1995, ApJ, 452, 156 
Cui, W. 2004, ApJ, 605, 662 

de Marco, D., Blasi, P., & Olinto, A. V. 2003, Astroparticle Phys., 20, 53 
Drury, L. O 'C. 1983, Rep. Prog. Phys., 46, 973 
Fossati, G., et al. 2000, ApJ, 541, 153 
Fraschetti, F., k Melia, F. 2008, MNRAS, 391, 1100 
Gabuzda, D. C. 1999, NewAstron. Rev. 43, 691 
Gaisser, T. K., Halzen, F., & Stanev, T. 1995, Phys. Rep., 258, 173; [Erratum: 1995, Phys. Rep., 271, 

355] 
Giovannini, G., Feretti, L., Venturi, T., Cotton, W. D., & Lara, L. 1999, in ASP Conf. Ser. 159, ed. 

L. O. Takalo & A. S. Silanpaa (San Fransisco ASP), BL Lac Phenomenon, a conference held 22-26 

June 1998 in Turku, Finland, 439 
Giroletti, M., Giovanni, G., Cotton, W. D., Taylor, G. B., Perez-Torres, M. A., Chiaberge, M., & 

Edwards, P. G., 2008, A&A, 488, 905 
Halzen, F., & Hooper, D. 2002, Rep. Prog. Phys. 65, 1025 
Halzen, F., Landsman, H., & Montaruli, T. 2005, arXiv: astro-ph/0503348 
Halzen, F. 2005, Lectures presented at the International WE Heraeus Summer School on Physics with 

Cosmic Accelerators [arXiv: astro-ph/0506248] 
Halzen, F., & Hooper, D. 2005, Astroparticle Phys., 23, 537 
Halzen, F., & Zas, E. 1997, ApJ, 488, 669 
Honda, M., & Honda, Y. S. 2002, ApJ, 569, L39 

15 



Honda, Y. S., & Honda, M. 2004a, ApJ, 613, L25 

Honda, M., & Honda, Y. S. 2004b, ApJ, 617, L37 

Honda, M. 2004, Phys. Rev. E, 69, 016401 

Honda, M., & Honda, Y. S. 2005, ApJ, 633, 733 

Honda, M., & Honda, Y. S. 2007, ApJ, 654, 885 

Honda, M. 2008, ApJ, 675, L61 

Honda, M. 2009, ApJ, 706, 1517 

IceCube Collab. 2004, Astropart. Phys., 20, 507 

loka, K., Razzaque, S., Kobayashi, S., &: Meszaros, P. 2005, ApJ, 633, 1013 

Kataoka, J., et al. 2001, arXiv: astro-ph/0105029 

Kazimura, Y., Sakai, J. I., Neubert, T., & Bulanov, S. V. 1998, ApJ, 498, L183 

Kotani, T., Kawai, N., Matsuoka, M., & Brinkmann, W. 1996, PAS J, 48, 619 

Koyama, K., Petre, R., Gotthelf, E. V., Hwang, U., Matsuura, M., Ozaki, M., & Holt, S. S. 1995, 

Nature, 378, 255 
Krawczynski, H., Coppi, P. S., Maccarone, T., & Aharonian, F. A. 2000, A&A, 353, 97 
Lagage, P. O., & Cesarsky, C. J. 1983, A&A, 118, 223 
Lagage, P. O., & Cesarsky, C. J. 1983, A&A, 125, 249 
Marscher, A. P. 1987, in Superluminal Radio Sources, ed. J. A. Zensus & T. J. Pearson (Cambridge: 

Cambridge University Press), 280 
Meisenheimer, K., Roser, H. -J., Hiltner, P. R., Yates, M. C, Longair, M. S., Chini, R., &: Perley, R. A. 

1989, A&A, 219, 63 
Miicke, A., & Protheroe, R. J. 2001, Astroparticle Phys. 15, 121 

Miicke, A., Protheroe, R. J., Engel, R., Rachen, J. P., &: Stanev, T. 2003, Astroparticle Phys. 18, 593 
Nan, R. D., Zhang, H. Y., Gabuzda, D. C, & Inoue, M. 1999, PASJ, 51, 955 
Nishikawa, K. -I., Hardee, P., Richardson, C, Preece, R., Sol. H., & Fishman, G. J. 2003, ApJ, 595, 

555 
Novak, G., et al. 2003, ApJ, 583, L83 

O'Sullivan, S., Reville, B., & Taylor, A. M. 2009, MNRAS, 400, 248 
Perley, R. A., Dreher, J. W., k Cowan, J. J. 1984, ApJ, 285, L35 
Plan, E., et al. 1998, ApJ, 492, L17 
Quinn, J., et al. 1996, ApJ, 456, L83 
Rawlings, S., & Saunders, R. 1991, Nature, 349, 138 
Rybicki, G. B., & Lightman, A. P. 1979, Radiative Procceses in Astrophysics (New York: Wiley- 

Interscience), ch.6 
Rebillot, P. F., et al. 2006, ApJ, 641, 740 
Sambruna, R. M., et al. 2000, ApJ, 538, 127 

Sembay, S., Edelson, R., Markowitz, A., Griffiths, R. G., & Turner, M. J. L. 2002, ApJ, 574, 634 
Silva, L. O., Fonseca, R. A., Tonge, J. W., Dawson, J. M., Mori, W. B., & Medvedev, M. V. 2003, 

ApJ, 596, L121 
Stecker, F. W., & Salamon, M. H. 1999, ApJ, 512, 521 
Takahashi, T., et al. 1996, ApJ, 470, L89 

16 



Tavecchio, F., &: Maraschi, L. 2001, in AIP Conf. Proc, 599, X-ray Astronomy: Stellar Endpoints, 
AGN, and the Diffuse X-ray Background, ed. N. E. White, G. Malaguti, & G. G. C. Palumbo, 
(New York: AIP), 979 

Toptygin, I. N., & Fleishman, G. D. 1987, Ap&SS, 132, 213 

Vietri, M. 1998, Phys. Rev. Lett., 80, 3690 

Waxman, E., k Bahcall, J. 1997, Phys. Rev. Lett., 78, 2292 

Woschnagg, K., for the AMANDA Collaboration, 2005, Nucl. Phys. B, Proc. Suppl. 143, 343 



17 



